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Abstract 

We consider the ensemble of real symmetric random matrices ifO.p) ob¬ 
tained from the determinant form of the Ihara zeta function of random graphs 
that have n vertices with the edge probability p/n. We prove that the nor¬ 
malized eigenvalue counting function of H^ n,p i weakly converges in average as 
n, p —> oo and p = o{n a ) for any a > 0 to a shift of the Wigner semi-circle dis¬ 
tribution. Our results support a conjecture that the large Erdos-Renyi random 
graphs satisfy in average the weak graph theory Riemann Hypothesis. 


1 Ihara zeta function of graphs and random matrices 

The Ihara zeta function (IZF) associated to a finite connected graph T = (V,E) is 
defined at u € C, for |u| sufficiently small, by 

ZrW = n (! - ^r 1 , (1.1) 

[C] 

where the product runs over the equivalence classes of primitive closed backtrack- 
less, tail-less cycles C = (ou, a. 2 , ■ ■ ■, on = aq) of positive length l in T, a* € V and 
z/(C) =1 — 1 is the number of edges in C [24], Being introduced by Y. Ihara [13] 
in the algebraic context, IZF represents now an intensively developing subject of 
combinatorial graph theory with applications in the number theory and the spectral 
theory (see e.g. [12, 24] and references therein); it has also been studied in various 
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other aspects, in particular in relations with the heat kernels on graphs [3], quantum 
walks on graphs [22], certain theoretical physics models [27]. 

Ihara’s theorem [13] proved first for (q + l)-regular graphs T says that the IZF 
(1.1) is the reciprocal of a polynomial and that for sufficiently small |u| 

Z v (u)~ l = (1 - u 2 ) r ~ 1 det(J + u 2 (B - I) - uA), (1.2) 


where A = (aij)i,j=i ) ...,n is the adjacency matrix of T = T n , n is the number of 
vertices of r n , B = diag(J2'j= i a ij)i=i,...,n and r — 1 = Tr(B — 2/)/2. This statement 
has been proven also for possibly irregular finite graphs (see [1, 23] for the combina¬ 
torial proofs of (1.2) and references related). The right-hand side of (1.2) represents 
an entire function; this means that Zp(u ) has a meromorphic continuation to the 
whole complex u-plane. Note also that r — 1 can be expressed in terms of the Euler 
characteristic of the graph because Tr (B — 2/)/2 = \E\ — n, where \E\ is the total 
number of edges of T. 

While the Ihara’s determinant formula (1.2) gives a powerful tool in the studies 
of the Ihara zeta function, the explicit form of Zr(u) can be computed for relatively 
narrow families of finite graphs. Regarding the case of infinite graphs, the definition 
of the Ihara zeta function represents an important problem that requires a number 
of additional restrictions and assumptions (see, in particular, [4, 10, 11, 17]). In 
the most cases, the graphs under consideration have a bounded vertex degree (in 
particular, regular or essentially regular graphs). 

A complementary approach is represented by a stochastic point of view, when 
the graphs are chosen at random from a set of all possible graphs on n vertices. This 
description ”in the whole” naturally leads to the limiting transition of infinitely in¬ 
creasing dimension of the graphs, n —>• oo. Certain aspects of large random d-regular 
graphs have been studied in [8]. The present note is related with the Ihara zeta func¬ 
tion (1.1) of random graphs whose average vertex degree p infinitely increases in the 
limit n — y oo. 

Let us consider an ensemble of n x n real symmetric matrices A^ n ' !p > whose entries 
are determined by a collection of jointly independent Bernoulli random variables 
An’ 1 = , 1 < i < j < nj such that 



1 — 5ij, with probability ^ , 

0, with probability 1 — ^ , 


(1.3) 


where 5ij is the Kronecker (5-symbol and 0 < p < n, p £ N. The adjacency matrices 
{R( n,p )} represent the ensemble of random graphs {r( n,p )} that can be referred to as 
the Erdos-Renyi random graphs [2], With this definition in hands, one can determine 
the random Ihara zeta function Z r ( n , P ) and to study it in the limit of infinite graph 
dimension, n —>• oo. In the present note, we consider the asymptotic regime of sparse 
random graphs, when !</)<«“ for any a > 0. 
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Passing to the normalized logarithm of (1.2), we get the following relation, 


log Z rM (u) = ^ Tr ( BM - 2/) log(l - u 2 ) 

logdet ((1 - u 2 )I + u 2 B {n 'ti - uA^ p ^ , (1.4) 

where 

n 

M.-fiiE"!?’' 1 

l] i=i 

Regarding the first term of the right-hand side of (1.4) 

0 (n,p)( u ) = _L Tl . ^ B (n, P ) _ 2 /) l 0 g(! _ 

it is easy to compute its mathematical expectation E0( n,/, l(u) with respect to the 
measure generated by the family An'- 1 (1.3), 

= (^—~p — 1^ ln(l — u 2 ) —•> (J^ — lj ln(l — u 2 ), n —»• oo. (1.5) 


The last expression shows that to obtain a finite value of E0( n, ^(«) in the limit 
n, p —>• oo, one has to rescale the parameter u = u p as follows, 


u 


2 

p 



Then the last term of (1.4) takes the form 


( 1 . 6 ) 


E^’ p \v) = i logdet ((1 - v 2 /p)I + H^P\v)^j , (1.7) 

where 

2 

H^p\v) = —B^ - —A^p\ ( 1 . 8 ) 

P Vp 

The presence of the factor p -1 / 2 in front of A^'p) is fairly natural from the point 
of view of the spectral theory of large random matrices. The normalization of bA>p) 
by v 2 /p is less common and arises because of the condition (1.5). Therefore it would 
be natural to say that {H^ n,p \v)} is the random matrix ensemble of the Ihara zeta 
function of random graphs. It is interesting to note that similar rescaling of the 
spectral parameter u 2 (1.6) by the vertex degree q + 1 is needed when instead of 
|p( n .p)| the ensemble of (q + l)-regular random graphs is considered in the limit of 
infinite q [18]. 
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Let us note that in the case of u = 1, the last term of the right-hand side of (1.4) 
gives the discrete version of the Laplace operator Ap 1 ’^ = — A^ n ' p ' 1 determined 

on the graph T (n >D . The spectral properties of Ap 1 ’^ of large random graphs have 
been studied in a number of works (see e. g. [5, 15]). The eigenvalue distribution of 
a version of the graph’s Laplacian given by p -1 / 2 Ap 1 '^ for finite and infinite values 
of 0 < p < n was studied in the limit n —> oo in [16]. Let us stress that this 
matrix p -1 / 2 Ap 1 ’^ is again different from (1.8). The ”zeta-function” normalization 
induced by (1.5) essentially changes the properties of the ensemble with respect to 
the graph’s Laplacian. In particular, in contrast to A r („, p ) the matrix H^ n,p ^ (v) with 
given v > 0 (1.8) is no more positively determined for large values of p, even the 
finite ones. 

Summing up, we can say that the random matrix ensemble H^ n,p \v) (1.8) is the 
new one that up to our knowledge has not been studied. In the present note we 
consider (1.8) with v € M and study the eigenvalue distribution of real symmetric 
random matrices H^ n,p \v) in the limit 

n, p —>• oo, p = o(n Q ), for any a > 0 (1.9) 

that we denote by (n, p) a —>• oo. Our main proposition is that the normalized 
eigenvalue counting function of H^ n,p \v) weakly converges in average to the well- 
known Wigner semicircle distribution shifted by v 2 . This statement is proved in 
Section 2. In Section 3, we discuss our results in relation with the limiting values of 
the correspondingly re-normalized Ihara zeta function for complex v € C. 

2 Limiting eigenvalue distribution of H^ n,p \v) 

Let us rewrite definition (1.8) in the form 

H^ p \v) = v 2 B^ - vA^ p \ »6l (2.1) 

Denoting the eigenvalues of H- n,rp (v) by Xi‘' p \v) < ■■■ < X^’ p \v), one introduces 
the normalized eigenvalue counting function, 

4"’ P) (A) = i# {j : \ ( "' p) (v) < A} , Ael. (2.2) 

The moments of this measure satisfy the following relation, 

M^’ p \v) = E QtV = J + X k dai n ’ p \X), k = 0,1,2,... , (2.3) 

where represents the averaged eigenvalue counting function, o4 r, ’ p) = Ecri n ’ p \ 

The main result of the present note is as follows. 
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Theorem 1 . For any given k € N, the averaged moment (2.3) converges in the 
limit (1.9), 


lim Mjf 1 ’^ (v) = fhk ( v ) 

(n,p) a —>-oo 



if k = 21 + 1, 


if k = 21, 


(2.4) 


where l = 0,1, 2 ,... and 


(2 p)\ 

lp p\(p + \)V 


P = 0,1,2,... 


(2.5) 


are the Catalan numbers. 

Proof. To study the moment , 


we consider the product 


4"’ p) (P, Q) = ^Y1 E (^ Pl " ' A Pa B qs 


2 Q = 1 


2020 


( 2 . 6 ) 


where Ylt=i Pi = -f*, 'f2t=i Qi = Qi P+Q = k, and where we denoted P = (p\,... ,p s ) 
and Q = (q ±,..., q s ). In (2.6), we assume p 2 > 1,... ,p s > 1, gi > 1, • • •, g s -i > 1 
and pi > 0 ,q s > 0. In what follows, we omit the bars over P and Q when no 
confusion can arise. 

We study the limiting behavior of variables L^ l,p \P,Q) with the help of the 
diagram technique close to that developed in [16]. We consider the product 

n 

(A Pl ---B q A = V (i Pl ). . (P (9l) ). . ■■■(A Ps ) / ^ . 

\ / 2020 * ^ ^ 202p-^ ' ' ' ' l p\-\ -H“l?20 ^ ' 2020 

io,ii,-..i Ps =l 


as a sum over the generalized trajectories (lp x \ Jql^ ,..., Tp} , Jqf ), where the closed 
trajectory of P steps is given by 

Pp i Pp2 ’ * * * ’ P { pJ ) (fOi ll •> • ■ • ; Ipi i l'p\ +1; • • • i 'Ipi P’P2 l • • • l^P— 1 5 ^C)) (2* 7) 

and = (j[ k) ,... k = l,...,s. 

Regarding (2.7), we associate to io a root vertex a and draw the new vertex 
f3, 7 , 5,... each time when we see a value of ii that is not equal to one of the values 
previously seen. We draw the blue edges that correspond to the steps of Pp and get 
a graphical representation of the trajectory Ip by a closed chain Gp with P blue 
edges. Clearly, Gp is a connected multi-graph with the root vertex a. We denote 
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by /3\,...,f3h the vertices of Gp that correspond to different values of variables 

*o, ipi > *pi+p2 5 • ■ • > ^pi+p2H—kps—1 1 1 — h < P■ 

The gi-plet Jq]^ can be represented by a set of q\ oriented red edges (/?i, j[) with 
1 < l < Qi- We can put a flesh at the head vertex 7 j. If the value of does not 
coincide with any element of Zp (2.7), then we say that the corresponding vertex 
7 r is the red one. In the opposite case we have € V(Gp). Then we construct, 
by the same procedure, a representation of the remaining parts Jq^\ ■ ■ . ,Jqf* and 
get a multi-graph that we denote by P(P, Q) = Gp hi f q- We say that Gp and 
represent the blue and the red sub-graphs of the diagram ZL(P , Q), respectively. 

Let denote by G = Gp and F = Fq the simple graphs associated with G = Gp 
and F = F Q . respectively and consider a part Fq of the red sub-graph of H(P,Q) 
such that there is no edge of E(F'q) that coincide with the elements of E(Gp). We 

_ V _ _ 

denote by V r {F') the set of red vertices with fleshes of Fq and write that Fq = 
Fq \ Fq . It is clear that 


n,,B(ip.J o ) = 0 IE(a)l+|B(f,)l 


where F' = Fq. Also it is easy to see that 


\C(G p ^Fq)\ = n(n- 1) • • • (■ n-\V(G)\-\V r (F')\+l ) = n IV(G)|+|l4(F')l(l+o(!)), n 00 . 


Then 


1 

np p /2+Q 


y. n a,b(z p ,j q ) 

{I p ,Jq}£C(G p VF q ) 


p\E{G)\—P/2 p \E(F')\-\E(F')\-\E(F")\ 

= n \E(G)\-\V(G)\+l n \E{F')\-\V r (F')\ ^ + n ^ °° ! ^ 2 ' 8 ) 

where F" = F'- y 

By using (2.8), it is not hard to prove that the terms of (2.6) that do not vanish 
in the limit (n, p) a —b 00 (1.9) have the diagrams that satisfy inequality 


\E(G)\ - \V(G)\ + \E(F')\ - \V r (F')\ + 1 < 0. (2.9) 

Indeed, if (2.9) is not satisfied, then the right-hand side of (2.8) gets a factor n~ k , 
k > 1 that suppress any power of p. The following two conditions are also necessary 
to have a non-zero limit: \E(G)\ > P/2 and 


\E(F')\ - \E(F')\ - \E{F")\ > 0. (2.10) 

Let us denote a diagram that verifies these three conditions by "H(P, Q ) = Gp tfcl Fq. 
Due to the Euler relation for the planar embedding of connected graphs, the only 
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equality is possible in (2.9) and this means that the simple graph G l+l F' is given 
by a tree. This tree is a plane rooted tree such that its blue sub-graph G is also a 
tree with P/2 = p edges. Indeed, let us denote by (<4, ei) the first leaf of the tree G. 
Then in the trajectory Ip (2.7) there is a couple of steps (ii-i, ii), (ii,ii+ 1 ) such that 
ii—i = ii+i- Removing this couple, we get the reduced trajectory X' p _ 2 such that its 
representation G is again a tree. Proceeding by recurrence, we see that P is pair, 
P = 2p and that \E(G)\ < P/2. Thus, p = \E[G)\. 

Let us note that given a graph Gp that is a rooted tree T pi one can easily 
restore the original multi-graph Tp and corresponding sequence of vertices W- 2 P by 
considering the chronological run over the tree T p . In this case the walk Wjjp is such 
that Tp is a multi-graph, where each couple of vertices { 7 , 5} is joined by either zero 
edges or exactly two edges - in there and back directions, ( 7 ,5) and (5,7). 

Regarding the red sub-graph F = F' U F" of 7~L(P,Q), one can easily see that 
\E(F')\ - \E(F')\ G {0, -1, —2,... } and that _\E(F")\ G {0,1, 2,... }. Then (2.10) 
is possible only when F" is empty and when F' = F’. The last relation means that 
the red part F' = F has no multiple edges. 

It is known that the number of rooted trees of p edges is given by the Catalan 
number t p . The positions and numbers of red edges being determined by (P, Q ), we 
deduce from ( 2 . 8 ) that the following relation is true, 

L ( //' P \P, Q) = t p v 2Q+p (l + o(l)), (n, p) a —>• 00 . (2.11) 

Returning to the moments (2.3) and regarding the trace 

1 i n 

- E (M,t) =-£ 

* 0 — 1 k times 

we have to choose 2 p elements H n of the last product that will be represented by 
—A. This can be done in ( I ways and this choice determines uniquely the s-plets 

\2p) 

P, Q. Combining this observation with (2.11), we get relation (2.4). Theorem 1 
is proved. □ 

It is known that the Catalan numbers t k verify the recurrence 

k 

4 +i = k > 0 and t 0 = 1 , ( 2 . 12 ) 

j =o 

and that the family of moments 

v 2p t p = f X 2p dp v (X), p = 0 , 1 , 2 ,... 

Jr 
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uniquely determines an even measure fi v with the density 


duv _ f f w _ 1 [V4v 2 - A 2 , if A G [~2v,2v], 

d\ 2irv 2 |o, otherwise, 

known in the spectral theory of random matrices as the semi-circle or the Wigner 
distribution [26]. It follows from (2.4) that the limiting moments rhk(v) can be 
represented as follows, 

rh k (v) = [ (v 2 + X) k d/j, v (X), k = 0,1,2,... (2.14) 

Jr 

and therefore the corresponding measure a v such that rhk(v) = / X k da v (X) is given 
by a shift of the semi-circle distribution (2.13), 

<r'(A) = - (A - v 2 ) 2 , \X — v 2 \ < 2\v\. (2.15) 

It follows from (2.15) that 

rhk(v) < (2v + v 2 ) k , k > 0 . 


The family of moments {rhk(v)} satisfies the Carleman condition and therefore the 
measure a v is uniquely determined. Thus, Theorem 1 implies the weak convergence 
in average of measures to a v \ this means that for any continuous bounded 

function / : 1 h>R the following is true, 

lim f f (X) dai n ' p \X) = [ f{X)da v (X). (2.16) 

n,P~yoo j R j R 

In this connection it should be noted that the generating function 


°° i p2v 

/(0=_ S 


dfly (^) 

—2v ^ £ 


verifies the following equation that can be easily deduced from ( 2 . 12 ), 

/«) = 1 


(2.17) 


-£-v 2 f{0' 

Regarding the generation function g(£) of rhj~{v) and using (2.14), we get equality 


1 ~ i \ f °° da v (X) 

9(0= "S{ 


k=0 


dg v ( X — v 2 

A-e 


= f(H-v 2 ). 



(2.18) 


This means that g{f) verifies the deformed version of (2.17), 

s(0 — 2 7 2 777' 

Relation (2.18) shows that the numbers rhk(v),k > 0 are determined by the recur¬ 
rence (cf. (2.12)) 

fc-i 

m fc+ i(u) = v 2 fh k (v ) + u 2 ^ m fe _i_.,-(u) rbj(u), fc > 1 (2.19) 

l=o 

with the initial conditions roo(v) = 1 and mi(o) = u 2 . In the particular case of 
u 2 = 1 the first ten values of nik(l),k > 0 (2.19) are given by 1, 1, 2, 4, 9, 21, 51, 
127, 323, 835. One can compute this also with the help of (2.4). These numbers are 
the Motzkin numbers [20]. 

Using the diagram technique developed above, one can get the following im¬ 
provement of Theorem 1. 

Theorem 2. Given k, the following relation holds, 

lim p - m k (v)) = R^\v), (2.20) 

(n,p) a —>oo V / 

where 



( 2 . 21 ) 

Proof. Let us consider a diagram 'H(P,Q) such that corresponding sum over 
C{'H{P, Q )) represents the term of the order 0{p~ l ) of in the limit n,p —>• oo, 

p = o(n). It follows from (2.8) that the diagram T~L{P,Q) verifies condition (2.9) 
and therefore its blue part G together with its red part F' represents a plane rooted 
tree G tfcl F' = In these relations, we denoted by G and F simple graphs obtained 
from corresponding multi-graphs G and F. The red part F' contains such red edges 
that do not coincide with edges of G. 

The next consequence of (2.8) is that the following condition is verified by hi, 

\E(G)\ - P/2 + \E(F')\ - \E{F')\- \E(F")\ = -1, (2.22) 

where F" = F \ F'. It is clear that (2.22) is possible in one of the three following 
situations: 
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- either \E{G)\ - P/2 = -1, \E(F')\ - \E(F')\ = 0 and \E{F")\ = 0, 

- or \E(G)\ - P/2 = 0, \E(F')\ - \E(F')\ = 0 and \E(F")\ = -1, 

- or |£(G)| - P/2 = 0, \E(F')\ - |£(F')| = -1 and \E(F")\ = 0. 

In the first case the chain of vertices of Gp is such that the corresponding walk 
W' 2 p generates a multi-graph of the tree-type such that there exists one edge passed 
four times and that remaining p — 2 edges are passed two times in there and back 
directions. The number of such walks is given by the formula 


# {m P } 


( 2 p)'- = p(p - 1 ) 

(p- 2 )!(p + 2 )! p (p + 2 ) 


(2.23) 


(see [14] for the proof). Then the remaining k — 2p red edges are to be distributed 
over 2p + 1 instants of time. This gives the factor ( 9 ^) in the right-hand side of 
( 2 . 21 ). 

The second case describes a graph such that the chain of vertices such that the 
blue diagram Gp is a rooted tree of p = P/2 edges passed two times and there exists 
one red edge of the form (/?, 7) such that 7 € V(Gp). It is not hard to see that 
obtain the corresponding chain, we have to consider a tree T p , to choose an edge e 
from it, to join one red edge 1 / to one or another side of e, to choose the orientation 
of the flesh of e! and to distribute k — 2p — 1 red edges over 2p + 2 instants of time. 
This gives the second term of the right-hand side of (2.21). 

Finally, the third case determines such diagrams T~L that Gp is a tree of p edges 
and there exist two red edges that make a multi-edge and the remaining red edges 
are simple. Let us describe how to construct corresponding diagrams and chains of 
vertices (walks). We consider two red edges of the form e' = (/5, 7 ) and e" = (/3, 7 ). 
Then we attach to the vertex f3 a tree 7 a that have a edges and point out an instant 
of time that will represent the starting and ending point of the corresponding walk. 
This can be done in 2a + 1 ways. Finally, we attach to /3 a tree 7 / of b edges, 
b = p — a in the way that it will represent the part of the walk performed between 
the fist and the second passages of (/?, 7 ) given by e' and e", respectively. Now it 
remains to distribute r — 2 — 2p red edges among 2p + 1 + 2 instant of times. This 
gives the factor ( 2 ^+ 2 ) • Taking into account elementary equality 


y] (2a + l)t a t(, — (p + l)tp_|_i, (2.24) 

a,b>0,a+b=p 


we get the third term of the right-hand side of (2.21). Relation (2.24) can be proved 
with the help of the generating function /(£) [14]. Also one can observe that the 
right-hand side of (2.24) represents the number of Catalan trees of p+ 1 edges with 
one marked edge; this gives the right-hand side of (2.24). 

Theorem 2 is proved. □ 
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3 Applications to IZF 


Variable (1.7) can be rewritten with the help of (2.2) in the following form, 

E^ p \v) = - logdet (l p + v 2 B^ - vA^ p A = / log (l p + \)da^' p \\), (3.1) 

where we denoted \ p = 1 — v 2 /p. Then the convergence of IZF for a sequence 
of graphs can be reduced to the question of the convergence of the corresponding 
spectral measures. This approach has been used for the first time in the studies of 
IZF of infinite regular graphs in paper [11], 

It is known that the normalized adjacency matrix A^ n,p ^ = p - 1 / 2 A^ n,p ^ of the 
Erdos-Renyi random graphs has all eigenvalues, excepting the maximal one, concen¬ 
trated with probability 1 in the limit n, p —>• oo on the interval [— 2 , 2 ], while this 
maximal eigenvalue is of the order ^fp [ 6 , 9]. Regarding (3.1) for the negative values 
of v, one can therefore expect that the limit 

lim E^ p \v) = E{v) (3.2) 

(n,p)a;—>-oo 


exists with probability 1 for all v such that —1/2 < v < 0. Moreover, observing 
that R is asymptotically close to the unit matrix in the limit (1.9), it would be 
natural to assume that convergence in average 

lim -E(log Z r( „, p) (v/y/p)) = ^~-E(v), (3.3) 

{n,p)a —>oo 71 Z 

where 

1 f 2 ^ , - 

E(v) = - — 5 - / log(l + A) \J 4u 2 - (A - v 2 ) 2 dX (3.4) 

2tTV 2 J—2\v\ 

holds for all v € (—1,1). While (3.3) is close to the weak convergence result 
—> a v (2.16) established by Theorem 1 above, one cannot use it directly be¬ 
cause the function log(l + A) is not bounded in the vicinity of —1. To justify (3.3), 
the detailed analysis of fine properties of eigenvalues of j s needed. This goes 

beyond the scope of the present note. 

Expression (3.4) can be rewritten in the form 


2 r 1 _ 

E(v) = — j log (l + v 2 + 2 vv) \/l — v 2 du, 


(3.5) 


where the last integral converges for any real v. Moreover, the right-hand side of 
(3.5) can be continued to a function holomorphic in any domain 


C € = {v : v € C, |u| < 1 — e}, 0 < e < 1. 
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Thus, assuming that (3.2) and (3.3) hold and remembering the relation between 
parameters u and v ( 1 . 6 ), one can say that our results support the conjecture that 
the family of random graphs satisfies in the limit (n, p) a —>• oo (1.9) a version 

of the weak graph theory Riemann Hypothesis [12, 25] that says that the Ihara zeta 
function Zy(u) is pole free for |it| < 1 /y/q, where u € C and q + 1 is the maximum 
degree of all the vertices of T. In our setting, the maximum degree normalization 
is naturally replaced by the averaged vertex degree p. Let us say that to prove this 
conjecture, one has to establish convergence (3.2) for complex v and this problem 
goes also far beyond the frameworks of our results. 

In the present note we follow a stochastic approach to the studies of the Ihara 
zeta function of graphs. We assume the graph’s edges to be present at random and 
consider the ensemble of such graphs in the limit of their infinite dimensions. Such a 
reasoning is completely in the spirit of the spectral theory of large random matrices 
used for the first time in the spectral theory of heavy atomic nuclei [ 21 , 26], see 
also [19]. In the frameworks of the stochastic approach, the collective properties of 
complex systems remain valid while the difficulties are mostly related with special 
cases that are relatively rare. This can be regarded as a kind of simplified description 
that nevertheless catches the principal features of the system under consideration. 

It should be noted that even in the frameworks of the stochastic approach to 
IZF of graphs, the problem of the limiting transition of the infinitely increasing 
dimension of graphs, n —>• oo, still persists. The main difficulty in the establishing 
of convergence of the normalized logarithm of IZF (1.4) is that because of possible 
presence of negative eigenvalues of one can guarantee the existence of the 

function zS n ’ p \v) (1.7) for small values of |u| only, and this smallness can converge to 
0 as n —>• oo. However, the proportion of the graphs that exhibit such anomalously 
small area of convergence can be sufficiently weak. The fact that the integral of the 
right-hand side of (3.5) converges confirms this conjecture. We also see that due to 
(3.3), a kind of the averaged version of the the graph theory Riemann Hypothesis 
can be true for the large irregular random graphs. 

Let us finally point out that statements similar to our results could be obtained 
for the family of d-regular random graphs . Regarding log Zp( n d) (u) (1.4) with 
the spectral parameter u = v/\Jd — 1 , we get the following version of (3.1), 

£ (n ’ d \v ) = — log det ( 7(1 + v 2 ) - ; V i (n ’ d A , (3.6) 

n V Vd- 1 / 

where A^ n ' d ^ is the adjacency matrix of With the help of the results of 

[18], one can show that the right-hand side of (3.6) converges as n, d —>• oo for all 
— 1 < v < 0 to the corresponding integral over the shifted semicircle distribution 
(2.15). Here again, regarding convergence of z^ n,d \v) in average, one can extend the 
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above domain up to v € (—1,1). To justify this, one would need to know the fine 
properties of the graph’s spectrum, including statements similar to the well-known 
Alon’s second eigenvalue conjecture (see e. g. [7]) as well as estimates of the lowest 
eigenvalue of A ( ' n,d ^j\Jd — 1 in the limit n, d —>• oo. 
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